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ABSTRACT 

Beraus** of <lic ir.ary co*.pUx subsystem Inter- 
actions vililch occur within nodem electronic poa^r- 
proccsslng systens. It Is often quite difficuu to 
accurately predict total system perfcraiance. A 
digital currpute** sImuUtlon tcduilque, wtilch ran 
be used to study such coi^poslte po/er-coftoltlonlng 
systCMS, Is applied to a spacccrcft dlrect-energy- 
transfer pen/er- processing systtsii. Tne results ob- 
tained duplicate actual system prrfrrmance with 
considerable eccuracy. The validity of the ap- 
proach and Its usefulness In studying various 
aspects of systr.i ^•rfo^m^fltc such as steady-state 
characteristics and transient responses to severely 
varying operating conditions are dcr.tOnst rated 
experimentally. 

INTRODUCTION 

During tne past two decades, the field of 
power electronics has been evolving and expand- 
ing Into an Increasingly i.ore Important element 
of modern electro- technology. New and i.iore power- 
ful methods of analysis and design have been dn- 
velO(>ed, resulting In Improved power-condition- 
ing equipment which often finds nrw and diverse 
applications In many areas of eU*ct»*or1cs. Re- 
cently, for reasons of economic efficiency and 
system reliability, there has been considerable 
Interest in developing standardized power-condi- 
tioning circuits and systems for space applica- 
tions. These standardized designs would obviate 
the need for co: tly, but often still not optimum, 
custom designs fer every new power system re- 
quired []]. Soiie svandaidizcd po./cr-con dltion- 
W. i-'dulcs have already been developed which 
s. 1,*^ a wide range of mission requirf.iicnts, 
thus I llmlfiating much of this cost and at the 
same time making available to the system d.signcr 
a pool of fully tested equipment with proven ro- 
llablllty [2,3]. 

trying to meet specific mission require- 
jinents with such standardized modules, It becomes 
Increasingly more Important to be able to study 
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the performance of power systems as a i/hole, 
rather than the various subsystems Individual Iv. 
Because of the many subsystem Interactions within 
such composite sysUns, It Is oft^n quite dif- 
ficult to predict the total system response even 
though the behavior of the Individual suhclrciiits 
may be familiar and well oocumenled. Mo.^ever, the 
behavior of the pwor-condl tioninn sysUm ar a 
whol^ when in steady-state operation or w:ien sub- 
jected to transient disturbances such as step 
changes In load or variations In the primary and 
secondary source characteristics must be fully 
Investigated before the Integrity of the aggregate 
power-condi tioning system caii be ensured. It is 
the purpose of this paper to present one approach 
to this problem of studying aggregate por#er-cof»d1- 
tlonlng systems, and to demonstrate the validity 
and usefulness of this approach. 

In considering the size and complexity of 
modern power-conditioning systcn»s, the digital 
camputer quickly comes to mind as an attractive 
and viable tool for such a study. Even tho«inh 
much progress has been made recently in developing 
and applying pov/erful analytical tools to the study 
of particular converter and Inverter circuits [4,5J, 
none of these techniques Is applicable to a general 
composite power-conditioning system which can be oVs- 
crlbed by an r.th-order system of nonlinear differ- 
ential equations. On the other hand, the rapid de- 
velopment of digital computer circuit and system 
analysis programs over the past few years has 
given the olectronics engineer the ahilUy to ex* 
amine more precisely very large analysis problems 
which had previously been Intractable or loo 
costly to pursue. Particular computer studies, 
both uigital and analog, have also been pre- 
sented recently which enhance U»e understanding 
of and eld In the design of individual ;)owcr con- 
verter circuits [6,/, 8], but again there have 
been no published results concerning simulations 
and studies of entire pov/er-condilloning systems. 
This paper describes In considerable detail the 
application of a gem ral **d1qital analog simula- 
tor" program, derived from IBM's Continuous System 
Modeling Program (CSMP) [9], tc e study of an 

aggregate pi»*ver-condtticning s> m being de- 

veloped at the NASA Goddard Space Might Center 
' for use on the International Ultraviolet Explorer 
(lUL) Spacecraft. 
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F1g. 1. Examples of the types of functional "blocks" 
.available for lirplomentlnn system descriptions In 
the computer program format. (A) Typo Hi mode- . 
controlled Integrator. (B) Type C, clipper. 

(C) Type n, relay. (D) Type +, weighted suimier. 


THE SIHULATIOIl TOOL 

As mentioned above, the tool used In this 
work Is an Interactive digital computer program 
based on an IBM application program, CSHP, as 
developed for use on the IBM 1130 system. The 
computer studies presented in tills paper were 
run on a Digital Equipment Corporation PDP-11/15 
digital computer equipped with a Tektronix <1013 
storage-tube graphics console using a substantially 
■ modified version of the 1130 CSHP. With appropri- 
ate modifications to the input data, however, sim- 
ilar studies could be conducted on a number of 
other computing machines 'ncluding the IBH 360/370 
systems. 

The. structure of this program is such that 
the particular dynamic system being investigated 
is entered into the computer in terms of tho 
system state equations. These equations are im- 
plemented by means of an appropriate connection 
of input/output functional relationships, much 
as is done when using an analog computer. Thirty 
types o-p functional "blocks" are available for 
this purpose, four of which are illustrated in 
Fig, 1. These symbolic blocks display all of 
the Information required to completely program 
a system simulation. Each block is given a 
single-character type designation, e.g. H, C, 

R or +, as shown in Fig. 1, which is uniquely 
defined within the program and corresponds to 
a particular function. Each block also must be 
assigned an individual idantilTcation number n. 
Tlie numbering of the blocks need not be sequen- 
tial since the program uses a sorting algorithm 
to determine the proper order for evaluating the 
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Fig. 2. Example of the application of the Simula 
tion tool to a second-order nonlinear system. 
(A) Van der Pol's equation. CB) Block diagram 
for solution of equation. Cc) Program input 
statements for output shown in (D). (D) Output 

data in form of X-Y plot. (H) Output data in 
form of Y vs time plot. 


output of each block. Specifying the input con- 
nections to each block of the system represen- 
tation and the associated parameter values com- 
pletes the simulation diagram. Depending on the 
block type, from zero to three inputs and from 
zero to three parameters must be specified. The 
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Fin. 3. Block diaprom representation of the DET 
system under investlpatlon with bold lines in- 
dlcatinp power channels. 

connections to the Inputs, shown syrboUcally In 
Fig, 1 by input variables e], en and C 3 , cor- 
respond to outputs from numbered blocks within 
the simulation diagram. The parameters, symbol- 
ized by PI, VZ and P3, correspond to real constant 
numbers which are a part of tl\e function repre- 
sented by the block. 

To clarify this procedure, a simple example , 
of a system represented by a second-order non- 
linear differential equation is shovm in Fig. 2. 

In this example, blocks 2 and 3 are standard in- 
tegrators, designated type I, which integrate the 
input signal continuously from the initial time 
to the final time, The initial condition for 
each Integrator can be specified as parameter PI 
of the respective blocks. In this example, the 
initial conditions have been specified as 0.0 and 
1.0 for blocks 2 and 3, respectively, as shown in 
Fig, 2(C). Blocks 4 and 5 are multipliers, desig- 
nated type X, which simply multiply one input by 
the other. The parameters PI, P2 and P3 associ- 
ated vHth block 1, are given values of -1,0, -c 
and c respectively, where c is a numerical con- 
stant of the system. For the case shown in Fig. 
2(C), c has been given a value of 5, All of the 
input data described above are entered into the 
computer in fixed format which greatly facilitates 
this routine task. 

After selecting an appropriate integration 
interval and a final solution time, the computer 
begins integrating tlie programmed system state 
equations using a second order Rungc-Kutta inte- 
gration scheme. The numerical value of the output 
of each of the blocks comprising the system repre- 
sentation is available as a tabulated printout 
and/or a graphical display at each integration in- 
terval. The speed with which the computer can 
perform these integrations is a function of the 
system complexity and the size of the integration 
step specified. Systens with short time constants 
may require relatively small Integration steps In 
order to accurately compute fast transients, and 
therefore would be subject to relatively long run 
times if the desired system response extends for 
a long period of time. Some experimentation with 
the integration interval is usually required to 
achieve a satisfactory combination of accuracy 
and computational speed. 



Fig. 4.‘ Decision flow chart illustrating the oper 
ational modes possible within the DET system, 
and the conditions v/hich activate them. 


Two of the most attractive features of this 
program are its interactive capabilities as orig- 
inally designed for the IBM 1130, and the flexi- . 
bility provided it, the choice of the output display 
formats as designed for the PDP-11/45 computer 
and Tektronix 4013 graphics unit. The interactive 
nature of the program enables the user to inter- 
rupt a simulation run at any time to modify the 
system description and thus perform on-line exper- 
imentation, Six'^een switches on the computer con- 
trol panel ore used for this purpose whereby such 
items as system parameters, initial conditions, 
and run controls can be adjusted. Additionally, 
the user may choose to display the output data of 
a given run in the form of an X-Y plot, with as 
many as five system variables plotted versus 
any other system variable with time an im- 
plicit parameter; or he may choose to display as 
many as five of the system variables as functions 
of time. In either case, numerical printouts of 
the same five or five different output variables 
are also available. Tv/o sample output displays 
for the example of Fig. 2 are presented in Fig. 

2(D) and (E). 

THE SYSTEM UNDER IHVESTIGATIOH 
The system being simulated with this program 
is the main bus regulator for the International 
Ultraviolet Explorer Spacecraft which is a joint 
mission hotvicen the Goddard Space Flight Center 
(GSFC/HASA), the British Science Research Council, 
and the European .Space Research Organization, 

This system is a direct energy transfer (DET) sys- 
tem whereby the primary source of electrical energy 
a solar array in this case, is coupled through a 
main distribution bus directly to the spacecraft 
electrical loads. The various power-conditioning 
subsystems, shown in Fig. 3, are activated only as 
needed, thus requiring the system to process only 
that amount of power which is needed to maintain 
the power bus at the specified voltage levei. The 
flov/ chart in Fig. 4 indicates the various modes 
of operation which are possible within this system. 
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Ftg. 5. Control unit outputs SI, to ,ost con- 
vertor, and S2, to charpo/shunt driver, as a 
function of bus voltaqc. 

The specified main bus voltage range for this 
mission is 28,»0.5G volts. The pm/er system control 
unit continuously monitors this voltage evel and 
generates appropriate signals to activa’i one of 
five different inodes of operation dcsigr,ed to meet 
this regulation requirement under widel;’ varying 
operating conditions. This control unic function 
is depicted graphically in Fig. 5. For a voltage 
of 28+0.14 volts, the system is said to bo in the 
deadband mode whereby all of the electrical energy 
required by the spacecraft loads is available from 
tlie solar array with no additional energy from 
the array to be dissipated and no supplementary 
energy from the secondary source required. Under 
these conditions, signals Si and S 2 from the 
control unit lie below the one-volt deadband 
threshold level, as shown in Fig. 5, and all of 
the subsystems remain idle. If, hov/ever, the 
solar array is unable to adequately supply the 
loads, the bus voltage falls and the control 
unit generates signal S] which activates the 
booa; mode of operation, and drives the boost 
converter with an adaptive duty-cycle control 
signal as determined by the secondary source 
voltage and the amount of additional power re- 
quired. The secondary source of energy in 
this case is a pair of rechargeable nickle- 
cadmiiim storage batteries which have terminal 
voltages between 17 and 25 volts, depending 
on the state of charge. If, on the other 
hand, the solar array is generating energy at a 
faster rate than is required by the spacecraft 
and the bus voltage rises above 28.14 volts, the 
excess energy must be Tunneled off to prevent an 
over-voltage condition on the main bus. Under 
those conditions, the control unit, in conjunc- 
tion with the charge/shiint driver and the battery 
charger subsystem, channels an appropriate amount 
of this energy into the secondary storage battery. 
Analog signals corresponding to the battery state 
of charge and the amount of excess power avail able 
from the solar array are combined to provide the 
duty-cycle control of the battery charger which 
ensures the proper charging rate. This mode of 
operation is referred to as the ehaAge. mode. If 
still more power is available than is required 
to properly cliarge the battery, the bus voltage 
continues to rise and additional energy must be 
dissipated, thus activating the menc 

of operation. When the bus voltage exceeds 28.18 
volts, the shunt elements are activated and dis- 



Flg. 6, Dissipative shunt regulator current as 
a function of bus voltage. 

sipate the excess power through resistors which 
are mounted on the spacecraft structure. The 
array of shunt elrnicnts is capable of dissipating 
up to a total of 450 watts of power as a linear 
function of the bus voltage as shown in Fig. 6. 

When the solar array is generating urmoaded energy, 
but at the same time tiie battery is fully charged, 
the fifth mode of operation is acti',et«d. In this 
case, all of the unneedod energy ift dissipated in 
the shunt pesistors and the system is said to be 
in the i/iuni medc. 

On examining the schematic diagram of this 
system as shown in Fig. 7, ont quickly recognizes 
that classical circuit analysis techniques may 
not bo satisfactory tools for studying the com- 
plex interactiens of the various subsystems, 
nevertheless, a careful examination of switching 
transients which n.ay occur on the main bus as 
a result of severely varying operating conditions 
is a necessary part of verifying the compliance 
of the system to its performance specifications. 

Of particular concern in a system such as this 
is its response to uncontrolUblo external inputs 
such as the intensity of sunlight on the solar 
panels or the state of charge of the secon- 
dary storage battery. Thus, the technique 
described in this paper for simulating the en- 
tire power system os a whole and systematically 
computing transient responses and steady-state 
characteristic*-, is a very jiiofi table approach 
and enables thr system engineer to more thor- 
oughly investigate his system vdiile it is still 
in the design stage. 

SYSTEM MODEL AND MATHEMATICAL I’.EPRESEHTATION 
The Model 

In order to effectively study the physical 
system described above, it is necessary to develop 
an appropriate model which is simple enough to 
yield a tractable mathematical representation but 
at the same tine detailed enough to enable an ec- 
curate description of those aspects of the system 
behavior which are of interest. It is clear that 
any given system c-an bo modeled in many different 
ways do'+'^nding on which system properties are con- 
sidered .>:ore important or of roost interest. Fre- 
quently, however, it is difficult to immediatoly 
determine the full significance of each component 
of the system, and therefore it often is necessary 
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F1g, 8. Idealized model of the power’-handllnq 
subnotwot'h of the boost converter subsystem. 
Actual component labels from Fig. 7 are shown 
next to the correspondlnri component model. 


to revise the system model at various stages of Its 
development when new evidence Indicates that In- 
correct assumptions have been made or that unneces- 
sary details which unduly complicate the model have 
been Included. By such a process, an appropriate 
model of the direct energy transfer system described 
In the preceding section v/as developed. The complete 
model consists of a tentli-order nonlinear network 
which approximates the power-handling subnetworks 
of the system, and several resistive Input/output 
relationships which represent the system control 
functions. 

As an example of the modeling approach taken 
In this v;ork, consider the schematic diagram of 
the boost converter subsystem as shoi/n In Fig. 7. 

All of the pov;er-handl1ng subnetworks of this sys- 
tem are drawn with bold lines at the top of the 
diagram. The boost converter subsystem consists 
of all of the components labeled v<ith a 500- series 
nunber. The power subnotivork Is shov/n in the up- 
per left-hand portion of the diagram, while the 
control subsystem Is shovm in the lower third of 
the figure. The model developed for the power- 
handling portion of this subsystem Is shown In 
Fig. 8. Each of the capacitors has been modeled 
as an Ideal lumped capacitance In scries v/lth a 
small Ideal lumped resistance which approximates 
the equivalent series resistance (ESR) of the 
actual component. Similarly, small winding re- 
sistances have been Included in the models of 
the Inductors and the transformer, but the winding 
capacitances have been neglected. The autotrans- 
former has bean modeled as a linear dgvico defined 
by the flux density versus magnetic field strength 
characteristic shown In Fig. 9{A). The 5em1con-_ 
ductor diodes have been approximated with a static 
nonlinear resistive characteristic as defined in 
Fig. 9(B), Similarly, the transistors are con- 
sidered to be cither fully on or fully off 
for the purposes of this study, and each Is mod- 
eled as an open circuit for the off condition 
and as a small saturation resistance for the on 
condition as Illustrated In l-'lg. 9(C). It should 



Fig. 9. Idealized models for (A) magnetic core of 
transformer T, (8) semiconductor diodes H3, D^-» 
and 05, and (C) transistor 02 and Q3 of Fin, 8. 

be noted that If detailed Infortmtlon concerning 
the performance of the semiconductor devices In 
this system Is desired, the model of each semi- 
conductor component can be changed to Include 
those properties of the device which «ru of par- 
ticular interest. The "ork described 1n this 
paper, however, Is con erned primarily vHth the 
transient behavior of the main bus voltage, and 
the dynamic performance of the semiconductor cci> 
ponents has not been considered. The voltage 
sources Vi and Vg shown In Fig. 0 represent the 
voltages induced across control v/indings 9-10 
and 12-11 respectively on transformer T504 which 
Is shown at the top of Fig. 7, These induced 
voltages arc assutiad tu constant during the 
on-time of the appre < 1' transistor and are 
therefore modeled as '.mt voltage sources. 
Resistor f ?5 Is a smal * ... I'cent sensing resistor 
which provides a current limiting signal to the 
control unit of this subsystem. 

' The control unit for the boost converter 
subsystem is modeled as ? purely resistive net- 
work. The purpase of this subnetwork is to pro- 
vide the duty-cycle control of the two pov/er 
transistors shewn In Fig, 8. This Is accomplished 
by generating a constant-frequency triangular f;ave- 
fora with a fixed peak-to-peak amplitude and add- 
ing to It a threshold level, Vj, v/b1ch 1s a func- 
tion of the bus voltage as shown graphically In 
Fig, 10. As the bus voltage falls, tlie threshold 
level falls, allff-ilng the power transistors to 
turn on for a greater portion of each cycle of the 
triangular wavefori. This Increased duty cycle 
eoobles the converter output current to incrocre 
to a new level which satisfies the Increased lead 
requirement. 

Af can be seen In Fig. 7, the various 
control circuits of the complete power system 
are often quite conplex, and the transfer func- 
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Hons derived ore by no menns 1r.ined1ately ob- 
vious. In most cnscs o substudy was conducted 
to reveal the salient features of these control 
subcircuits and l ) precisely .dotermine the appropri- 
ate transfer characterisitic. Some of these 
studies involved experitientnl measurements, 
while in other cases digital computer circuit 
analyses were run to reveal the significant 
circuit relationships, In all cases, hov;evcr, 
the assumption tiiat these circuits are essentially 
purely resistive vtas carefully verified, 

in a similar manner, models have been developed 
for the battery ctiarger subsystem, the shunt elements, 
the charge/shunt driver, and the main control unit. 

The complete system model includes controllers to 
doteraiine the dufy-cycle of the pov/or transistors in 
the boost converter and the battery charger subsys- 
tems, and a controller to determine the value of the 
shunt dump resistance at all times. 

The Hathcniatical R e presentation 

To analyze the model presented in the preced- 
ing subsection, it is necessary to construct a 
inatharatical representation for it. As described 
above, the control subcircuits of this system have 
been modeled as purely resistive networks. Thus, 
the behavior Of these subnetoorks is completely 
specified by a set of nonlinear functional equa- 
tions which may be derived from the laws of the 
elements in these circuits and the laws of inter- 
connection of them} or they may be approximated 
from experimentally measured transfer characteris- 
tics of the circuits. Given these equations, all 
of ‘he voltages and currents in the control sub- 
networks can be computed under all possible operat- 
ing conditions. But in practice, only those com- 
binations of curriiits and voltages vdiich c-xplicitly 
determine the control signals of interest have been 
represented mathematically. The coefficients of 
these functional equations are determined from 
appropriate combinations of the actual circuit 
component values, including variable elements such 
as the battery state of charge, bus voltage, and 
solar array output current. Thus, the currents 
and voltages of interest in these circuits, and 
therefore the system control signals, change value 
as pertinent system variables change; and these 
changes in control signals give rise to the vari- 
ous modes of opP**ation of the system as described 
in the third i,ection of this paper. 

To complete the mathematical representation 


of the madol, ib'* poi/''r-lmndlinq subortworks must 
bo considi red, iitue the toriiincd suhryst'i i have 
been nci.'.*1i*d os a itntli-ordi'r nonllrit'ur dyriiinlc r.i*!- 
v/ork, j, set of ten ititli?|'i'n'l(;nt nonliii^Mr differentitil 
equations is required to iMthematically specify 
the behavior of the nnlv.wk. Having developed 
these stale equations from the laws of the ele- 
ments and the lows of interconnection of the net- 
work, and given the values of the state variables 
at initial time tn, the currents and voltages 
tlirougtiout this nclv/orl: con bo detcnmlnod for any 
time t > I q. But because of the many different 
operating conditions pnssible within this system, 
one set of state cquitions is not sufficient to 
completely characterize the system at all titnes. 
These different modes of operation are the result 
of 0 transistor or a diode switching on or off , 
and thus effectively changing the topology of the 
pov/er-handling subnetworks and consequently the 
laws of interconnection of them. Therefore, a 
sot of state equations must ho written for each 
possible mode of operation within the system, and 
the state variables must be continuous through 
a transition from one mode to another os dictated 
by the control functions; i.e., the final value 
of state ^ for one mode of operation rust be the 
Initial state for the next mode. With these 
sets of differential equations and the set of 
functional equations described above, the entire 
system model is accurately represented mthemati- 
cally at all times and for all operating condi- 
tions. The next section describes how this com- 
prehensive maUiomatical representation is Imple- 
mented in the digital computer. 

PROGRAIMIHG THE SYSTEM 

As discussed earlier in the paper, the program 
used in this work may be considered to be a digi- 
tal analog simulator in that the system solution 
is computed on a digital conputer, but at the same 
time the program input format resembles analog com- 
puter connection diagrams. Thus, given the roathe- 
matical representation of the system as presented 
above, an analog- type connection diagram can be 
drawn as an intermediate step to programming the 
system. An example of this procedure for one of 
the system state equations is shovm in Fig. 11. 

In Fig, 11(A) are the differential equations v/hich 
determine the trajectory of the state variable vci 
of the boost converter power subsystem through all 
possible inodes of operation of that subsystem. 

The connection diagram in Fig. 11(B) illustrates 
the i/ay in which each of those differential equa- 
tions is constructed from the pertinent system 
state variables and their associated coefficients 
v/hich are determined nunterically from appropriate 
combinations of the circuit component values. 

For exoniplo, the derivative of the voltage across 
capacitor Cl during the time transistor Q2 of 
Fig. 8 is on is a vieightcd sum of the currents 
through the inductors Lq and L|,, and the flux in 
transformer T. Block ni of Fig. 11(B) performs 
this summation where • 


PI « 

P2 = - 


1 

cr 
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ID) 

Flu. n. DlagroiRflnn the state equations, 

(A) Diffprential equations describing the tro- 

. 8 for three 
Block dia- 
gram) reprcscntatloh of these equations. 

The parameters i, u. and A represent, respectively, 
tlic Clean magnetic patli length, the penneablllty, 
and the cross-sectional area of the magnetic core 
of transformer T, Similarly, blocks n^ and 

compute the derivative of this st.itc variable for 
other operating conditions. Functional relay 
blocks are used to select the proper combination 
of variables at each Instant of time as deter- 
mined by the appropriate control function, thus 
simulating the switching semiconductor compo- 
nents of the system. In this case, control sig- 
nal Y determines If either transistor Q2 or Q3 
Is on 5 and 1f so, control signal X determines 
which one. These decisions are made at each 
computational time Increment, and thus the cor- 
rect value of the derivative of vci with respect 
to time Is available for all time and all opera- 
ting conditions. The output of relay block ng 
Is Integrated numerically to yield the continu- 
ous trajectory of the state variable vci. Con- 
tinuity from one isode of operation to another 
Is maintained through the memory of the Inte- 
grator; 1.0. t the last output of the Integrator 
for one operating condition becojnes the initial 
condition for the next mode. In a similar man- 
ner, the functional representations of the con- 
trol subsystems which generate the signals X 
and Y are diagrammed and connected to the ap- 
propriate relay blocks In the state equation 
diagrams. 

When the entire system Is d1.ngrairred In 
this manner, It Is then a simple matter to 
complete the prograriimlng procedure. After 
labeling each of the blocks in the connection 
diagram, the entire system description con be 
entered Into the computer through preformatted 
statements as Illustrated in Fig, 2(C). When 
the last block description Is entered, the In- 
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teractlve program asks tho user to enter the 
Integration Interval, the final solution tine, 
onT the output data opt1<ns and form-at. After 
processing these entries, the computer begins 
calculating the system state trajectory and 
displaying the output d5ta requested. The user 
can observe this d.ita as it Is being computed 
and has the option of interrupting tli" run at 
any time t*> tiwdify the system descrliitlon or 
output requests as suits his needs. In this 
v;ay, on-line cxpcrlr^ntation can be conducted 
as Is described In the follov/lno sections. 

EXPtRlKUITAL VERIFICATIOll 

To verify the ability of the computer repre- 
sentation described above to simulate actual sys- 
tem porfonnance, a series of system tests were 
conducted on the esscirblcd DET prototype unit and 
duplicated on tlie ccr.piftcr model. A comparison of 
the wosured and cojputod results from two of these 
tests Is presented In this section. The tests v/erc 
designed to activate each of the system modes of 
operation, and thereby enable an examination and 
evaluation of all phases of the computer model. 

The data recorded during tliese tests ore displayed 
In Figs. 12 through 17. The (A) portion of each 
figure Is an oscillogran of the experimentally mea- 
sured response of the actual OCT power-condltlcnlng 
system, while the (D) portion of each figure Is the 
corresponding graphical computer output for the sim- 
ulated system response. The coiiputer generated out- 
put has been photographically reproduced to yield 
a grid of the same vndth ns the exijerimi'ntal oscil- 
lograms, This moans that the horizontal time axis 
In parts (A) and (B) of the same figure are Identical. 
In each of these cscillcgrams and computer displays, 
the saiiKi four system variables— load current, shunt 
current, battery current and bus voltage--are plot- 
ted versus time with the same vertical scale fac- 
tors of volts or cmpercs per division. Those vari- 
ables and their associated seal big are presented 
In Table 1 for cosy refemince when examining the 
data displays. It should be noted that the coi putev* 
generated grid docs not have the same vertical 
to horizontal ratio as the oscilloscope grid; i.c., 
the actual spacing betvjcen major vertical divisions 
In part (B) of each figure Is greater than that In 
part (A), Thus, the various computer generated 
plots a; ; tar elongated In the vortical direction. 

The third v/avefori' In each of the displays is the 
secondary-battery current which has been desig- 
nated os positive when charging and negative when 
discharging the battery. 

The first test, prcfscntcd in Figs. 12, 13 
and 14, subjects the OET system to transient 
disturbances which cause it to go from on Initial 
steady-state condition 1n the charge/shunt mode 
of operation to the boost mode and back to the 
charge/shunt .node In a period of 10 milliseconds. 

This transient response Is caused by a step change 
In the system load current from 0 to 8 A at 1.7 
milliseconds and tack to zero at 6.7 milliseconds, 
while the solar array supplies 4 A of current 
continuously througheut the test. When no load 
current Is being tJrav.-n, the entire 4 A of solar 
array current must be shunted from the bus to 
avoid an over-voltage condition. During this 
test, the signals corresponding to the battery 
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state of choroe have been sot to indicate o discharged 
condition, and tlwreforo the battery is charged at a 
iiaxliniini rate of opproulmately 1.? A. The rewining 
Z.3 A of solar array current is shunted into the 
shunt resistors. When the system loads require 8 A 
of current, however, the boost converter must be 
activated to supply the additional 4 A of load current 
required. The complete transition from charfie/shtmt 
to boost and back to charge/shunt operation is illus- 
trated in Fig. 12, The horironal scale for this 
figure is 1 millisecond per major divisioni The 
oscillogratiB are riultiplo exposures, and therefore 
the switching of the load current for each of the 
four traces does not nocess. rily occur at precisely 
the same position on the grid. Also, there is some- 
times a zero offset in the measured results due to 
a drift in the oscilloscope between the multiple 
exposures. 

The first 1,7 milliseconds of Fig. 12 illus- 
tra' ' steady-state operation in tlie chargc/shimt 
me The bus voltage Is steady at a level of ap- 
proximately 28.2 volts. At 1.7 milliseconds, the 
required load current is switched from 0 to 8 A. 

This sudden change in load current causes a down- 
ward jump in the bus voltage due to the voltage 
drop across the equivalent series resistance of the 
output capacitor, rr? of Fig. 0, of the boost con- 
verter. This capacitor had been charged to the 
steady-state value of bus voltage, 28,2 V, and at 
the instant of sv/itching is the only source of 
current to supply the new load requirement, Uith 
this rapid drop in bus voltage to a new level cor- 
responding to the boost mode or operation, the 
shunt, elements and the battery charger turn off 
consecutively, and the boost converter turns on. 

This transient response is illustrated with an 
expanded time scale in Fig. 13, v/here the hori- 
zontal scale has been changed to 100 micro- 
seconds per maior division, and the switching 
occurs at the ,70 microseconds position of the 
trace. In this expanded view, it becomes clear 
that there is a delay of approximately 40 micro- 
seconds between the time the signal to turn off 
the shunt olecients is generated and the time 
that the shunt current actually falls to zero. 

This delay, duo primarily to the time required 
for operational amplifiers in the shunt drive 
circuitry to coma out of saturation, causes 
the bus voltage to continue to fall even as 
the boost converter supplies additional cur- 
rent to the loads. 

Looking again at Fig. 12, one sees that 
the shunt current does fall to zero and remains 
zero while the battery current decreases to a 
steady-state value of approximately -6 A, and 
the bus voltage settles to an average value 
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Fig. 12, Test Case I. Solar array current “ 4A. 
Load current sviitches from 0 to 0 to 0 amperes. 
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Fig. 13, Test Casa I, Expanded view of 0 to 8 
ampere svHtching transient. 
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Fig. 14. Tost Case 1. Expandod view of 0 to 0 
ampere switching transient. 


of approximately 27.0 V. At approximately 
6.7 milliseconds, the load current is switched 
back to zero, again causing a jump in the bus 
vc'tage, this time in the upward direction, duo to 
the equivalent series resistance of the boost convert- 
or output capacitor. This switching transient is 
expanded in Fig, 14 where again the expanded time 
scale is 100 microseconds prr major division. The 
bus voltage is seen to jump to a new level of approxi- 
mately 20.2 V v)hich demands that the boost converter 
turn off and that the energy stored In the output 
inductor of the boost convertor. Lb in Fig. 8, as 
well as the 4 A solar array current be shunted from 
the bus. 

As can be seen in Fig. 12, during the interval 
that the system is returning to a steady-state condi- 
tion in the charge/shunt mode, the battery current 
does not change from -6 A to 1.2 A in a simple singlo- 
time-constant manner. The shape of this transient 
vjaveform is a very complex nonlinear function of the 
state variables of the systeiii. As diodes and transis- 
tors are switched on and off and as energy stored in 
inductors and capacitors is depleted or increased, 
currents and voltages throughout the system undergo 
abrupt transitions and proceed along complex tfajcc- 
tories. It is beyond the intended scope of this 
paper to analyze this dynamic behavior in detail, 
but it should be pointed cut that the computer- 
simulation waveforms presented in this section dupli- 
cate with considerable accuracy the shapes of the 
actual system responses and consequently tin' complex 
system interactions which cause them. It is there- 
fore possible to use the computer simulation to 
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Fig. 15. Tost Case II. Solar array current OA. 
Load current sv/itches from 0 to 8 to 0 amperes. 


examine any of the state variables which are used 
in describing the system, including the transformer 
flux, to determine more precisely the v;ay in v/hich 
the various subsystems interact to produce these 
observed v/aveforms. With the computer simulation, 
it is a simple natter to look at any or all of the 
system state variable trajectories, whereas when 
working with a breadboard unit it ofter is difficult, 
and sometimes Impossible, to make accurate moas- 
ments of the variables of Interest. 

A second test case is presented in Figs. 15, 

16 and 17. Again, the same four vmefonns are dis- 
played in each figure with the vertical scales 
listed in Table 1, Throughout this test, the solar 
array is generating no current. The transient dis- 
turbance is again introduced by a step change in 
load from zero to 0 A and back to zero over a period 
of 10 milliseconds. Fig. 15 displays the entire 10 
millisecond response. For the first 1.7 and the final 
3.3 milliseconds of this display, no external load 
current is supplied. However, at all times the vari- 
ous control subsystems of this DET system t'oquire 
some power to operate continuously In a standby 
mode so as to be able to respond immediately to any 
system disturbance, l/hen no solar array current 
is available, this standby pov/cr must be supplied by 
the battery through the boost converter, thus 
t'cquiring the system to operate in the boost made. 

The amount of current- required to maintain this 
standby condition is approximately 10(1 inA and 
thus is not discernible in the graphic displays 
presented. The boost mode of operation is in 
evidence, however, when examining the bus voltage 
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Fig. 16. Test Case II. Expanded viev/ of 0 to 0 Fig. 17. Test Case If. Expended viev; of fi to 0 

amporc switching transient. aicnore sv/ltching transient, 


waveform In Fig. 15. Over the last three mini- tion, thus giving a fairly complete represen- 

seconds of this trace, the bus voltage slowly tatlon of Its perfori'ince. The oscillograms 

decays unti l the boost mode of operation Is octiv- and corresponding computer generated vaveforn: 

aled and the standby power Is supplied through have been presented in a form w.t1ch hopefully 

the boost converter, enables the reader to make a careful comparison 

and draw his ovm conclusions as to the accuracy 

At 1.7 milliseconds, the load current 1s of the slmulntlon tC' ‘miquo described In this 

switched from zero to 8 A and the system niovos . paper. The simulation runs illu tratc.1 In Figs, 
further Into the boost mode of operation. This 12(B) and 15(C) reqire approxlr.' *. ' ,■ ’ , i!.ir of* . 

switching Interval Is Illustrated v;lth an ex- computation tine, while each of tb nil- 

panded time scale of 100 microseconds per major second runs requ<.-«s approximately " i.ites 

division In Fig. 16. At 6,6 milliseconds the to complete the deta display. - i r 'id pre- 
load current Is switched back to zero and the viously, the small computer used In ’ work 

transient of Fig. 17 results. When the 0 A v)os a P0P-11/A5 with 16K words of me: . y using 

load Is removed from the bus, the bus voltage a disk operating system. The approx w. ate cost 

jumps upward to a value greater than 28 V as per hour of computation, including purchase price 

a result of the boost convertor current suddenly amortization and raintalnonce, v;as ten dollars, 

flowing through the equivalent scries resistance 

of the' output capacitor, and the chargc/sliunt COUCLUSIOHS 

n»t!e of operation Is momentarily activated. 

The bus voltage quickly falls back Inside the A digital computer simulation technique has 

deadband region, however, and then sloivly decays been presented as a viable tool for the study of 

to the upper limit of the boost node of oper- aggregate power-co.nditlonlng systems. As discussed 

ntlon In order to again supply the system stand- earllci In the paper, a preliminary step in the de- 

by po'.vcr. velopr.ent of such a system simulation Is the postu- 

lation of a conceptual model. This model a*id Its 

The presentation of experimental and cor- subsequent mathematical representation ore of utnost 

responding simulatlotr data In this section has. Importance in the simulation process, for it is at 

because of space limitations, been restricted tl)1s level that the nature of the real sysien holng 

to two test cases. These particular cases simulated is captured. The computer program performs 

v/cre selected from a large pool of data which the repetitive computations necessary to generate 

was collected throughout this endeavor because system state tralcctorics, but the computations per- 

the switched load tests cause the DET system, formed are dictated by the differential and functional 

and likewise Its computer simulation model, to equations developed to represent real system perforn- 

functlon in each of the possible modes of opera- ance. The resultant computer simulations of actual 




systcri trajectories can tijy be os iic'uroto os the 
systenn model allows, and only tlmsc aspctlo of tin; 
syslrti) behavior bhkh ore Included in tlie sli ulotloii 
model con be expectrd to aiipear in the coMputed re- 
sults, 

The examples presented above ore not intended 
to be exhaustive. This poper presents the icoults 
of the development of a tool for use in studying 
power-proct*5s1n<i systems which arc too complex to 
bo handled v/lth closslcol analysis tools such os 
phase plane and lin* orlzation techniques. Hav- 
ing verified the v.Oidlty and ot curacy of this tech- 3. 
nlqiie, It is intendtjd that the comjmtor model devel- 
oped bo fully exorcised to aid in the final devel- 
opment of the DCT system described above and to con- 
tribute, from the beginning stages, to the development 
of future poner-ccmditioning systemi. One of the pro- 4. 
posed uses for this program is a paromotor optinii* 
zntion study through repeated simulation runs to dc- 
torainc which v.if.iua of selot-l'nl coi.iponrnts night 
improve some aspect of the system perfortnanco. This 
parameter modification can be accorplislied quite 
easily within the simulation program, whereas con- 
siderable difficulty can sometintas he encountorod 
when attempting such a study at the breadboard level. 

Other uses of such a program include the study of 
semiconductor component stresses wi thin the system, 
especially during severe switching transients such 
as illustrated above } and the study of system 
stability when in any of the modes of operation 
possible. 

The digital computar technique described in 
this paper has been applied quite successfully 
to a particular dircct-cnergy-transfer pov.'er- 
conditioninq system of considerable complexity 
which is still in the development stupes, The 
computer simulation of this system has been ex- 
ercised vigorously and extensively In an attempt 
to verify its accuracy under all operating con- 
ditions and to demonstrate its usefulness. It is 
felt that the test data presented in this paper are 
a good representation of this testing. Moreover, 
it' Is believed that the functional block diagram 
onnroach used in this simulation technique possesses 
significant advantages over the network topology 
approach of other commonly available electronic 
circuit analysis pronrams when used to study such 
composite nonlinear dynamic systems. 
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